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Nonlinear theory of solitary waves 
associated with longitudinal particle motion in lattices: 
Application to longitudinal grain oscillations 
in a dust crystal * 
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Ruhr-Universitat Bochum, D-44780 Bochum, Germany 
(Dated: submitted: 17 dec. 2003, accepted: 6 jan. 2003) 

The nonlinear aspects of longitudinal motion of interacting point masses in a lattice are revisited, 
with emphasis on the paradigm of charged dust grains in a dusty plasma (DP) crystal. Different 
types of localized excitations, predicted by nonlinear wave theories, are reviewed and conditions 
for their occurrence (and characteristics) in DP crystals are discussed. Making use of a general 
formulation, allowing for an arbitrary (e.g. the Debye electrostatic or else) analytic potential form 
0(r) and arbitrarily long site-to-site range of interactions, it is shown that dust-crystals support 
nonlinear kink-shaped localized excitations propagating at velocities above the characteristic DP 
lattice sound speed vq. Both compressive and rarefactive kink-type excitations are predicted, de- 
pending on the physical parameter values, which represent pulse- (shock-)like coherent structures 
for the dust grain relative displacement. Furthermore, the existence of breather-type localized oscil- 
lations, envelope-modulated wavepackets and shocks is established. The relation to previous results 
on atomic chains as well as to experimental results on strongly-coupled dust layers in gas discharge 
plasmas is discussed. 

PACS numbers: 52.27.Lw, 52.35.Fp, 52.25.Vy 



I. INTRODUCTION 



A wide variety of linear electrostatic waves are known 
to propagate in plasmas 0, 0] ■ It is now established that 
the inherent nonhnearity of electrostatic dispersive me- 
dia gives birth to remarkable new phenomena, in par- 
ticular related to the formation and stable propagation 
of long-lived nonlinear structures, when a balance be- 
tween nonhnearity and dispersion is possible jSj, iJj . Since 
about a decade ago, plasma wave theories have received 
a new boost after the prediction (and subsequent ex- 
perimental confirmation) of the existence of new oscil- 
latory modes, associated with charged dust-grain mo- 
tion in dust-contaminated plasmas, as well as the pos- 
sibility for an important modification of existing modes 
due to the presence of charged dust grains 0^0. A 
unique new feature associated to these dusty (or com- 
plex) plasmas (DP) is the existence of new strongly- 
coupled charged matter configurations, held responsible 
for a plethora of new phenomena e.g. phase transitions, 
crystallization, melting etc., and possibly even leading 
to the formation of dust-layers (DP crystals) when the 
inter-grain potential energy far exceeds the average dust 
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kinetic energy; a link has thus been established between 
plasma physics and solid state physics 0. These dust 
Bravais-type quasi-lattices, which are typically formed in 
the sheath region in low-temperature dusty plasma dis- 
charges, and remain suspended above the negative elec- 
trode due to a balance between the electric and grav- 
ity forces 13, 0, Hol , are known to support harmonic 
excitations (acoustic modes) in both longitudinal and 
transverse-shear (horizontal-plane) directions, as well as 
optical-mode-like oscillations in the vertical (off-plane) 
direction 12, .13, 15, 18, 17,ii[li. 

The longitudinal dust-lattice waves (LDLW) are rem- 
iniscent of waves ('phonons') propagating in atomic 
chains, which are long known to be dominated by non- 
linear phenomena, due to the intrinsic nonlinearities of 
inter- atomic interaction mechanisms and/or on-site sub- 
strate potentials 20, 21, 22, J3,,^]. These phenomena 
have been associated with a wealth of phenomena, e.g. 
dislocations in crystals, energy localization, charge and 
information transport in bio- molecules and DNA strands, 
coherent signal transmission in electric lines, optical pulse 
propagation and many more [2^ |2^, |23, IH, ll^i] . Even 
though certain well-known nonlinear mechanisms, e.g. 
shock formation, electrostatic pulse propagation and in- 
stabilities, have been thoroughly investigated in weakly- 
coupled (gas-like) dusty plasmas [s^. El] , the theoreti- 
cal investigation of the relevance of such phenomena with 
waves in DP crystals is still in a pre-mature stage; apart 
from the pioneering works of MelandsofT^. who first de- 
rived a Korteweg-DeVries (KdV) equation 'si'l associ- 
ated with longitudinal dust-lattice oscillations, Shukla 
[T9|. who predicted the formation of dust cavitons due 
to lattice dynamical coupling to surrounding ions, and 
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the investigation of related nonlinear amplitude modula- 
tion effects by Amin et al. js^l a little later, not much 
has been done in the direction of a systematic elucidation 
of the relevance of dust-lattice waves being described by 
the known model nonlinear wave equations. It should, 
however, be stressed that some recent attempts to trace 
the signature of nonlinearity in experiments |34l l35l Is^ 
have triggered an effort to interpret these results in terms 
of coherent structure propagation [s^ |3^ [s^, ■ es sen- 
tially along the physical ideas suggested in Ref. |l2l|. 

In this paper, we aim at reviewing the procedure em- 
ployed in the derivation of a nonlinear evolution equation 
for longitudinal dust grain motion in DP lattices, and dis- 
cussing the characteristics of the solutions. Emphasis is 
made on the methodology, in a quite exhaustive manner, 
in close relation with previous results on atomic chains, 
yet always focusing on the particular features of DP crys- 
tals; we will discuss, in particular: 

- the physical assumptions underlying the continuum 
approximation; 

- the choise of truncation scheme, when departing from 
the discrete lattice picture; 

- the long-range electrostatic interactions, differenti- 
ating DP crystals from ordinary classical atomic chains 
(spring models); 

- the physical relation between different solutions ob- 
tained. 

Some of the results presented here are closely related to 
well-known previous results, yet enriched with a new an- 
alytical set of coefficients allowing for any assumed range 
of site-to-site interactions and any analytical form of the 
interaction potential. The present study is, therefore, 
valid in both short and long- Debye length DP cases, and 
also aims at providing a general 'recipe' which allows one, 
for instance, to assume a modified (possibly non-Debye- 
type) potential form and obtain the corresponding set 
of formulae in a straightforward manner. In specific, we 
have in mind the modification of the inter-grain interac- 
tions due to ion flow in the sheath region surrounding 
the dust layer, which may even lead to the crystal being 
destabilized, according to recent studies from first prin- 
ciples 39, 40J. 

Most of the results presented here are general and ap- 
ply, in principle, to a sufficiently general class of chains 
of classical agents (point masses) coupled via arbitrary 
(and possibly long-range) interaction laws. Neverthe- 
less, our specific aim is to establish a first link between 
existing nonlinear theories and the description of lon- 
gitudinal dust-lattice oscillatory grain motion in a DP 
crystal. At a first step, our description cannot help be- 
ing 'academic', and somewhat abstract: an ideal one- 
dimensional DP crystal is considered, i.e. a single, unidi- 
mensional, infinite-sized, dust-layer of identical (in size, 
charge and mass) dust grains situated at spatially peri- 
odic sites (at equilibrium). Effects associated with crys- 
tal asymmetries, defects, dust charging, ion-drag, dust 
mass variation and multiple dust-layer coupling, are left 
for further consideration ^41] . Transverse (off-plane) mo- 



tion, in particular, will be addressed in a future work. 



II. THE MODEL 

A. Equation of motion 

Let us consider a layer of charged dust grains (mass M, 
charge Q, both assumed constant for simplicity) forming 
a Bravais lattice, of lattice constant tq. The Hamiltonian 
of such a chain reads 
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where r„ is the position vector of the n— th grain; 
Unmifnin) = Q 4>{x) is a binary interaction potential 
function related to the electrostatic potential 4'{x) around 
the m— th grain, and r„,„ = |r„ — Vm\ is the distance 
between the n— th and m— th grains. We shall limit our- 
selves to considering the longitudinal x) motion of the 
n— th dust grain, which obeys 



M 



dt 



dx-n 



(1) 

where E{x) = —d(t){x) / dx is the electric field; the usual 
ad hoc damping term is introduced in the left-hand-side 
(Ihs), involving the damping rate i/, to account for the 
dust grain collisions with neutrals. Note that a one- 
dimensional (ID) DP layer is considered here, but the 
generalization to a two-dimensional (2D) grid is straight- 
forward. At a first step, we have omitted the external 
force term F^xti often introduced to account for the ini- 
tial laser excitation and/or the parabolic confinement 
which ensures horizontal lattice equilibrium in experi- 
ments 115. The analogous formulae for non-electrostatic, 
e.g. spring-like coupling interactions are readily obtained 
upon some trivial modifications in the notation. 

The additive structure of the contribution of each site 
to the potential interaction force in the right-hand-side 
(rhs) of Eq. ^ allows us to express the electric field in 
(HJ as: 



E{x) 



^ m 



^ ^ 1 d''+V(r) 

^ ^ l'\ dr^'+^ 
1=1 i'=i 



r=lro 



[{dXn+l - 6x„y - {SXn - SXn^lY 



(2) 



where I denotes the degree of vicinity, i.e. 1 = 1 ac- 
counts for the nearest-neighbour interactions (NNI) and 
I > 2 accounts for distant- (second or farther) neighbour 
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interactions (DNI). The summation upper limit N nat- 
urally depends on the model and the interaction mecha- 
nism; even though N 'traditionally' equals either 1 or 2 in 
most studies of atomic chains, one should consider higher 
values for long-range-interactions e.g. Coulomb or De- 
bye (screened) electrostatic interactions (the latter case 
is addressed below, in detail). In the last step, we have 
Taylor-developed the interaction potential 0(r) around 
the equilibrium inter-grain distance Ivq = |n — m|ro (be- 
tween /— th order neighbours), viz. 



1 d''0(r) 



l'=0 



r—\n—m\ro 



where I' denotes the degree (power) of nonlinearity in- 
volved in each contribution: /' = 1 is the linear interac- 
tion term, V — 2 stands for the quadratic nonlinearity, 
and so forth. Obviously, 5xn — Xn — Xn^ denotes the 
displacement of the th grain from equilibrium, which 
now obeys 



M 



d'^iSXn) 



d{SXn] 



dt^ dt 

Q{ (P'^ro) {SXn+l + SXn-1 - 2SXn) 



N 



(f>"{lro) {SXn+i + 6Xn-l - 26Xn) 



1=2 



1 d/+Hr-) 



l'=2 



l'\ dr^'+^ 



r=ro 



[[5Xn+l - SXnY - {SXn - SXn-lY ] 



N oo 



1 d/+\r) 



1=2 l'=2 



l'\ dri'+^ 



r=lro 



{SXn+l-SXnY -{SXn-SXn-lY 



(3) 



We have distinguished the linear/nonlinear contributions 
of the first neighbors (lst/3rd lines) from the correspond- 
ing longer neighbor terms (2nd/4th lines, respectively). 

Keeping all upper summation limits at infinity, the last 
discrete difference equation is exactly equivalent to 
the complete equation However, the former needs to 
be truncated to a specific order in I, I' , depending on the 
desired level of sophistication, for reasons of tractability. 



B. Continuum approximation. 

We shall now adopt the standard continuum approxi- 
mation often employed in solid state physics 7], trying 
to be very systematic and keeping track of any inevitable 
term truncation. We will assume that only small dis- 



placement variations occur between neighboring sites, i.e. 

dXn±l = OXn ± Itq— + ^il-^O) 



I.e. 



5Xn 



+1 



^^" = ^"°a^ + 2('"") dx^ 



E 



m! 9a;™ 



(4) 



and 



6x^-Sx^.,^lro-~-{lro) ^ 



m! dx"^ 



(5) 



where the displacement Sx{t) is now expressed as a con- 
tinuous function u = u{x,t). 

Accordingly, the linear contributions (i.e. the first two 
lines) in (O now give 



N 



25Xn) 



oo N 



Q2m^ 



qe 



m=l 1=1 
2 



- (2m) 



(2m)! 
Y.cl^"{lro){lro)'- 



N 



MY. 



1=1 



dx'^^' 



C2n 



dx"^ 



m 



— M I C2 Uxx "t" C4 Uxxxx ~t- C(5 tlxxxxxx •• 



(6) 



where the subscript in Ux denotes differentiation with 
respect to x, i.e. Uxx = d'^u/dx^ and so on. We see 
that only even order derivatives contribute to the linear 
part; this is rather expected, since the model (for v = 
0) is conservative, whereas odd-order derivatives might 
introduce a dissipative effect, e.g. via a Bur ger s-like (~ 
Uxx) additional term in the KdV Eq. below |42l li^ |43| . 
The definition of the coefficients C2m (m ~ 1,2,...) is 
obvious; the first term reads 



^2 = §-oE^"(M'' 



N 



M 



2 2 



(7) 



1=1 
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which defines the characteristic second-order dispersion 
('sound') velocity vq [cf. Vp in (6) of Ref. [s^], related 
to the longitudinal oscillation eigenfrequency wq.l; also 

1 O ^ 

1=1 



2 g 



1=1 



(8) 



and so on. Notice that vf = Vq/12 for NNI, i.e. if (and 
only if) one stops the summation at Imax — N = I, 
like Eq. (26) in Ref. HI [and unlike Eq. (5) in Ref. 
|38l |. whose 2nd term in the rhs is rather not correct, for 
I 7^ 1 i.e. DNI]. See that the 'relative weight' of any given 
2m— th contribution as compared to the previous one is 
roughly (2m - 2)!/(2m)!, e.g. 4!/6! 1/30 for m = 3, 
which somehow justifies higher (than, say, w = 2) order 
contributions often neglected in the past; nevertheless, 
this argument should rigorously not be taken for granted, 
as a given function u{x,t) and/or potential 4>{x), may 
present higher numerical values of higher-order deriva- 
tives, balancing this numerical effect; clearly, any trun- 
cation in an infinite series inevitably implies loss of in- 
formation. 

We may now treat the quadratic nonlinearity contri- 
bution in (PI (the last two lines for = 2) in the same 
manner. Making use of Eqs. 10} and (0), and also of the 
identity — 6^ = (a + h){a — b), one obtains 



1 ^ 



1=1 



{SXn+l-SXnf-{SXn-SXn-lY 



oc oo 



EE 



^ (2m- l)!(2m')! 

= lm' = l ^ I \ J 



N 



CXD OO 

^^EE 



d 



2m— 1^, Q2ra' 



U O U 



dx2m-l g^2m' 
Q2ra-1^ Q2m' ^ 



Qx2m-1 Q^2m' 



m — 1 m' — l 
M I Cli U^x + Cl,2 Ux 



+ C2,1 UxxUxxx + •••j • (9) 

The definition of the coefficients Cm.m' {n^, m' = 1, 2, ...) 
is obvious: the first few terms read 



Ci,i^%rlY^rilro)l', 



N 



M 



(10) 



1=1 



which defines the first nonlinear contribution [eg. B in 
Eqs. (5) and (7) in Ref. [s^; we note that a factor 1/2 



and 1/M is missing therein, respectively], 

N 



1=1 



1!4! M 



C2,l 



2 ^~-5^r(M^^ 



3! 2! Af " ^ 

1=1 



and so on. 

The cubic nonlinearities in Q (the last two lines for 
I' = 3) may now be treated in the same manner. Making 
use of Eqs. (@J) and |SJ|, as well as of the identity: a'^—b^ = 
(a — fe)(a^ + ab + b"^) , one obtains 



1 ^ 
^^E'^"'(^^o) 



1=1 



{6Xn+l-SXn)^-i6Xn-SXn-iy 



^ OO OO oo 

4EE E 



i-(-ir +(-1) 

(2m)! m'\ m"\ 



m' -\-m" 



X 



m—1 m' — l ra" — 1 

^ d^'^u d"''u 



^<?!)""(/ro) (Zro)2"+"'+™" 



1=1 



Q^2m Qj.ni' Q^m" 



my: 



Q2m^ Qm'^ Qm'\ 



-'m.rn' ,m" 



m—1 



g^2m g^m' Q^m" 



M 



Cl,l,l {UxY Uxx + (Cl,l,2 + Ci_2,l) Ux {u^xY 



+ Cl,2,2 {Uxxf + 



(11) 



The definition of the coefficients Cm,m',m" {m,m' ,m" = 
1,2,...) is obvious; their form is immediately deduced 
upon inspection, e.g. 

Cl,l^ = \%4Yr'{lr,)l\ 



1=1 



Cl,l,2 



-ci,2a 



12 M 



N 

riy^r'{lr,)l\ (12) 



1=1 



and so forth. We note that the second term in Eq. 111(1 
cancels. 

Higher order nonlinearities in Eq. Q (the last two 
lines therein for > 4), related to fifth- (or higher-) order 
derivatives of the interaction potential 0, will deliberately 
be neglected in the following, since they are rather not 
likely to affect the dynamics of small grain displacements. 
It should be pointed out that, rigorously speaking, there 
is no a priori criterion of whether some truncation of the 
above infinite sums is preferable to another; some ad hoc 
truncation schemes, proposed in the past, should only 
be judged upon by careful numerical comparison of the 
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relevant contributions - e.g. in Eqs. ©, ©, lfTT|l above 
- and/or, finally, a comparison of the analytical results 
derived to experimental ones. 

Keeping the first few contributions in the above sums, 
one obtains the continuum analog of the discrete equation 
of motion 

il + v u-vl u^x = vj rl Uxxxx - Po Ux + qo {uxf u^x , 

(13) 

which is the final result of this section. Notice that 
UxUxx = {ul)x/2] also, {uxYuxx = (-ui^)a;/3. The co- 
efficients 



phenomenological; the damping effect may then be re- 
inserted in the analysis at any step further, by plainly 
adding a similar ad hoc term to the equation(s) model- 
ing the grain dynamics. It may be noted that damping 
comes out to be weak, in experiments ^35j , so one may in 
pinciple proceed by including dissipation effects a poste- 
riori, and then comparing theoretical or numerical results 
to experimental ones. 



An exactly computable case - the Debye 
ordering 



C2 : 



2 2 

^^1 ^0 



C4 , 



-Cl,l 



90 = ci,i,i , 



defined by Eqs. 0, ©, (CDl) and (O, respectively, 
should be evaluated for a given potential function 0, 
by truncating, if inevitable, all summations therein to a 
given order Imax- Note that, quite surprisingly, the infi- 
nite neighbour contributions may be exactly summed up, 
in the case of Debye (screened) electrostatic interactions, 
as we shall show below. Let us point out that Eq. l(T^ 
is general; the only assumption made is the continuum 
approximation. Also, should one prefer to improve the 
above truncation scheme, e.g. by including more nonlin- 
ear terms, one may readily go back to the above formulae 
and simply keep one or more extra term(s); in any case, 
one can find the exact form of all (retained and trun- 
cated) coefficients above. On the other hand, Eq. H13() 
generalizes the previous known results for monoatomic 
lattices in that it holds for an arbitrary degree of inter- 
site vicinity (range of interactions). 

Let us point out that the above definitions of the co- 
efficients in Eq. (|13|l are inspired by the Debye-Hiickel 
(Yukawa) potential form (whose odd/even derivatives are 
negative/positive; see below), in which case they are de- 
fined in such a way that all of v^, vf, po and qq take 
positive values. Nevertheless, keep in mind that the sign 
of these coefficients for a different potential function (p 
is, in principle, not prescribed; indeed, analytical and 
numerical studies of the nature of the inter-grain inter- 
actions from first principles suggest that the presence of 
ion flow, for instance, may result in a structural change 
in the form of (p, leading to lattice oscillation instabil- 
ity and presumably crystal melting jl^; see e.g. Refs. 
[23,0. However, our physical problem loses its mean- 
ing once this happens; therefore, we will assume, as a 
working hypothesis in the following, that C2 and C4 bear 
positive values (so that vq , vi are real) - as a requirement 
for the stability of the lattice - and that, in principle (yet 
not necessarily), the same holds for po and qq. 

We observe that, upon setting i/ = 0, qq — 0, rg ~ a 
and I = 1 (NNI), which imply that vl ~ v'^a? /12 and 
Po = -Qa?(j)"'{a)/M = ^{a)a^lM in Eq. one re- 

covers exactly Eq. (26) in Ref. [ij [also see the definition 
in Eq. (16) therein]; also cf. Ref. |33|. Equations (5) - 
(7) in Ref. are also recovered. 

In the following, we will drop the damping term [second 
term in the right-hand-side of Eq. (|13|) ]. which is purely 



Most interestingly, the summations (in /) in the above 
definitions of coefficients Cm,m',... above, converge and 
may exactly be computed in the Debye-Hiickel (Yukawa) 
potential case: (pDif) — Qe^''/^°/r, for any given num- 
ber N of neighboring site vicinity: iV 1 for the near- 
est neighbor interactions (NNI), N — 2 for the second- 
neighbor interactions (SNI) and even N equal to infinity, 
for an infinite chain. The details of the calculation are 
given in the Appendix, so only the final result will be 
given here, for later use in this text. Note the definition 
of the lattice parameter k — ro/Xi), to be extensively 
used in the following; in fact, k is roughly of the order of 
(or slightly above) unity in laboratory experiments. 

Truncating the summations at = 1 (NNI) , relations 
(0, ®, ^ and (Hai give 



L,0 



(.r"'))V(.^A|,) = 12(.r^))V(..^A!,), (14) 



(NNI) 12Q 



These relations coincide with the ones in previous studies 
for NNI 15, "H]. 

Truncating the summations at A^ = 2 (SNI), relations 
0, ©, Cnil and ^ give 

( {SNI).2 _ f 1 + K + K,y2 ^ ^+K+H ^ 



D 



(17) 



accompanied by an extended set of expressions for 
^(SJV/)2 („(5^'f))2/i2, now, unlike in the NNI case 
above), p^^^^^ and q\^^^^ (see in the Appendix for de- 
tails). 
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For higher Imax = N, even though the effect of adding 
more neighbors is cumulative, since all extra contribu- 
tions are positive, these diminish fast and converge, for 
infinite N, to a finite set of expressions, which can be 
calculated via the identities: X^z^i ~ a/(l — a) and 
Si^i "^V^ — ~ln{l ~ a) (for < a < 1); details can 
be found in the Appendix. This procedure is similar to 
the one proposed in Ref. [3 and later adopted in Refs. 
mHi. One obtains 



2Q2 1 



2 \ 2 



— csch — 



ln(l 



(18) 

for the characteristic oscillation frequency wl,o = 
vq/ (kXd)', the result for vq is obvious; cschx = 1/ sinhx. 
A numerical investigation shows that the numerical value 
of the frequency in the region near rg ~ A/j (i.e. k « 1) 
is thus increased by a factor of 1.5 or higher, roughly, 
compared to the NNI expression above (see Fig. and 
so does the characteristic second-order dispersion veloc- 
ity Vq = LJL,o ro = LOLfl A_D n (see Fig. A similar effect 
is witnessed for the characteristic velocity vi, related to 
the fourth-order dispersion 



M\d 96k V 2 



{f? + 2) coshn + 2 (k^ — 1 + k sinhn) 
(see Fig. and for the nonlinearity coefficients 



Po 



MXd Ke'^ - 1)3 
+ e'^^{K^ + 3k +6) 



6 + e''(K^ -3k- 12) 



6 



ln(l + sinh k — cosh k) > (20) 



and 



MX, 



1 



(e« - 1)4 



-12 



-e^'-|lK 



12k + 48) cosh K 

+ 2(k3 - 6k - 18) + 2(k2 - 6) sinh k] 

12 1 
ln(l + sinhK — coshK) > . (21 



Upon simple inspection of Figs. 2|and[Sl one deduces that 
go takes practically double the value of pq everywhere, 
and thus draws the conclusion that go should rather not 
be omitted in Eq. (O (cf. e.g. Refs. El US HI HI ) , 
for the case of the Debye potential. 



III. LINEAR OSCILLATIONS 



Let us first consider the linear regime in longitudinal 
grain oscillations. For the sake of rigor, one may revert to 
the discrete formula ^ and consider its linearized form 
by simply neglecting the two last (double) sums therein. 
Inserting the ansatz 5xn ~ expz(nfcro — wi), where u is 
the phonon frequency and k = 27r/A (respectively. A) is 
the wavenumber (wavelength), one immediately obtains 
the general dispersion relation 



- iv] 



AD ^ 

^ Y: ^"{Iro) sin- 



I kvQ 



1=1 



— 2^ (IkXd) sm . (22) 

1=1 



One may readily verify that the standard ID acoustic 
wave dispersion relation uj fa vq k is obtained in the small 
k (long wavelength) limit: check by setting sin(Z fcro/2) « 
Zfcro/2 (and recalling the general definition of wq above). 
Of course, taking this limit simply amounts to linearizing 
the continuum equation of motion derived above (and 
keeping the lowest contribution in k). As pointed out 
before (see e.g. Ref. 0|), one thus recovers the dust- 
acoustic wave dispersion relation obtained in the strong- 
coupling dusty plasma regime (upon defining the density 



Ud as 



which may nevertheless appear somehow 



^19^ heuristic in this ID model). 

Notice that the form of the dispersion relation, in prin- 
ciple, depends on the value of N. However, in the case 
of the Debye interactions, i.e. explicitly substituting 
(P4)d{x) /dx"^ into H22|l . one obtains 



oj {to + iiy) = 



N 



4Q^ 



-i. 2 2k + (^k)^ 
(/k)3 



sm 



I kvQ 



(23) 

A numerical investigation, e.g. for k = 1 (see Fig. 1^, 
suggests that the dispersion curve quickly sums up to a 
limit curve, even for not so high values of N (practically 
for = 2 already). The values of the frequency reduce 
with increasing k, as suggested by the exponential term. 
We see that the dispersion curve possesses a maximum 
at fc = Tr/ro — vr/ {nXjj) for any value of k and N. 

The dispersion curves of dust-lattice waves have been 
investigated by both experiments (see e.g. Refs. [i^li^'l 
and ab initio numerical simulations [4^ . In should never- 
theless be acknowledged that the results of these studies 
do not absolutely confirm the dispersion curves obtained 
above, which suggests, as poined out in Ref. li^, that 
one-dimensional crystal models may be inappropriate for 
real dust crystals. 
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IV. THE KORTEWEG - DE VRIES (KDV) 
EQUATION 

In order to take into account weak nonlinearities, a 
procedure which is often adopted at a first step consists 
in keeping only the first nonlinear contribution in Eq. 
H13|l (by cancelling the last term in the rhs, i.e. setting 
go = 0) and then considering excitations moving at a 
velocity close to the characteristic velocity vq. A Galilean 
variable transformation, viz. 



: (24) 



then provides the Korteweg - De Vries (KdV) Equa- 
tion 



0, 



(25) 



where a term Utt was assumed of higher-order and thus 
neglected. The coefficients are 



a = 



2vo 



,,2^2 



2vo 



sgnpo =po/bo| 



(26) 

We have introduced the parameter s (~ +1 / — 1), denot- 
ing the sign of po, which may change the form of the solu- 
tions (see below); as discussed above, it is equal to s = -f 1 
for the Debye-type interactions. It should be noted that 
this procedure is identical to the one initially adapted 
for dust-lattice-waves in and then followed in Ref. 
IMj Jiij-2SJ (for s — +1) as may readily be checked, yet 
the new aspect here lies in the generalized definitions of 
the physical parameters above. Also notice that positive- 
oriented (~ x) propagation was considered; adopting the 
above procedure in backward (~ —x) propagation is triv- 
ial, yet it should be carried out by re-iterating the ana- 
lytical procedure and not by plainly considering v — > —v. 
the KdV equation is not symmetric with respect to this 
transformation (also see that the velocity v appears un- 
der a square root in the formulae). 

As a mathematical entity, the KdV Equation has been 
extensively studied ^ 2L M, 50, 51, 52, 53, 54, 5^, 
so only necessary details will be summarized here. It 
is known to possess a rich variety of solutions, includ- 
ing periodic (non-harmonic) solutions (cnoidal waves, in- 
volving elliptic integrals) |5J]. For vanishing boundary 
conditions, Eq. (|25|l can be shown (see e.g. in Refs. 
[27I Issj ) to possess one- or more (N—) soliton localized 
solutions u'Ar(C,T) which bear all the well-known soli- 
ton properties: namely, they propagate at a constant 
profile, thanks to an exact balance between dispersive 
and nonlinear effects, and survive collisions between one 
another. The simplest (one-) soliton solution has the 
pulsed-shaped form 



wi{(,t) — —s wi^msech^ 



(C 



Co)/Lo 



(27) 



where xq is an arbitrary constant, denoting the initial 
soliton position, and v is the velocity of propagation; in 



principle, v may take any real value even though its range 
may be physically limited, as in our case, where v has 
been assumed close to vq] this constraint will be relaxed 
below. A qualitative result to be retained from the soliton 
solution in H27() is the velocity dependence of both soliton 
amplitude wi^m and width Lq, viz. 

= ^v/a = 6vvo/\po\ , 



L,^{ib/vf^^[2vlrl/{vvo)]'/'. 

We see that wi^^Lq — constant^ implying that nar- 
rower/wider solitons are taller/shorter and propagate 
faster /slower. These qualitative aspects of dust-lattice 
solitons have recently been confirmed by dust-crystal ex- 
periments fs^. Notice that the solutions of (|25|l satisfy 
an infinite set of conservation laws 3, 55] ; in particu- 
lar, the solitons wn carry a constant 'mass' M ~ / wdC, 
(which is negative for a negative pulse), 'momentum' 
P ^ J w'^dC, 'energy' P S{w'i/'^ + 'U'^)dC,, and so 
forth (integration is understood over the entire x— axis) 
[5^ I55I I. See that the forementioned amplitude-width 
dependence of the 1-soliton solution l|77|l is heuristically 
deduced from the soliton 'mass' conservation law (im- 
plying conservation of the surface under the bell-shaped 
curved in Fig. [TJ: taller excitations have to be thinner 
and vice versa. 

Inverting back to our initial reference frame, one ob- 
tains, for the spatial displacement variable u{x,t), the 
kink/antikink (for s = —1/ + 1) solitary wave form 



ui{x,t) 



tanh 



{x — vt — xo)/Li 



(28) 



which represents a localized region of compres- 
sion/rarefaction (for s ~ +1/ — 1), propagating to the 
positive direction of the x— axis (see Fig. [7||. The am- 
plitude Mi,m and the width Li of this shock excitation 
are 



Ml.' 



Li 



ro 



bol 



2vf 



[2 Vq {v - Vq)] 



1/2 



1/2 



bol Ar, 



imposing 'supersonic' propagation {v > vq) for stability, 
in agreement with experimental results in dust crystals 
[35j |. Notice that faster solitons will be narrower, and thus 
more probable to 'feel' the lattice discreteness, contrary 
to the continuum assumption above; therefore, one may 
impose the phenomenological criterion: L ro, amount- 
ing to the condition v/vq 3> 1 + 2v1/vq [« 1.17 for the 
Debye NNI case; see (|14f) above], in order for the above 
(continuum) solution to be sustained in the (discrete) 
chain. Nevertheless, supersonic wave stable propagation 
has been numerically verified at a wide range of velocity 
values in atomic chains [23 . 15^ . where Eqs. H13|) and H25|) 
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arise via a procedure similar to the one outlined above; 
also see Ref. for a recent experiment in crystalline 
solids. Finally, note that vq in real DP crystals bears 
values as low as a few tens of mm/sec [33.l35| |. 

Remarkably, Eq. 125[l is exa ctly solved via the Inverse 
Scattering Transform |5.ll IH^l IHnj . for any given initial 
condition 0), which is generally seen to break-up into 
a number of (say N, depending on w(C,0) [5^) solitons 
plus a tail of background oscillations. These considera- 
tions, including, in particular, the two-soliton solution W2 
of Eq. (|25|l , which represents two distinct humps moving 
at different velocities and colliding during propagation 
without changing shape, have been postulated to be of 
relevance in the interpretation of recent dusty plasma 
discharge experiments js^ [s^ . 

The wide reputation of the KdV Equation H25() is 
mostly due to the exhaustive knowledge of its analyti- 
cal properties 0, IHIEl El El El El , 

in addition to 

its omni-presence in a variety of physical contexts, not 
excluding the physics of ordinary (ideal, i.e. electron- 
ion) plasmas 3, '4'| and, more recently, dusty plasmas 
jS, ,.3Q]. However, in the above dusty-plasma-crystal con- 
text, it has been derived under specific assumptions (low 
discreteness and low nonlinearity effects; also, a propa- 
gation velocity v ~ vq) which may be questionable, in a 
real DP crystal. Even if the first one is virtually impos- 
sible to cope with, analytically, the latter ones may be 
somehow relaxed via a different approach, to be outlined 
below. 



V. HIGHER-ORDER KORTEWEG-DEVRIES 
(EKDV) EQUATIONS 

In order to derive a KdV equation from the continuum 
equation of motion 1)13(1 . we have neglected the coeffi- 
cient q, which is related to the cubic nonlinearity of the 
interaction potential. Nevertheless, a simple numerical 
investigation shows that this term is not small, and may, 
in certain cases, even dominate over the quadratic term p, 
as in the Debye potential case (see the discussion above). 
Therefore, one is tempted to find out how the dynamics 
is modified if this term is taken into account. 



not necessarily, in general), so that the nonlinearity co- 
efficients, i.e. —sa (for s = +1) and a, bear negative 
and positive (respectively) values in this (Yukawa crys- 
tal) case. 

The EKdV Eq. was thoroughly studied in a clas- 
sical series of papers by Wadati 21], who derived it for 
nonlinear lattices, then obtained its travelling-wave and, 
separately, periodic (cnoidal wave) solutions and, finally, 
exhaustively studied its mathematical properties. Both 
compressional and rarefactive solitons (say, W2,±-, to be 
distinguished from the KdV solution wi) were found to 
solve Eq. (for either signs of s); adapted to our 

notation here jS^], they are of the form 



'.v/{Ccosh^[{C-vT-Co)/Lo] 

+ Dsmh'[iC~VT~Co)/Lo]}, (30) 



and 



(Cr) = +sv/{Dcosh^[iC ~VT- Co)/io] 

-fCsinh2[(C-«r-Co)/io]}, (31) 



where 



C 



\2v^ 

a 
6 



6av 



- 1 



= ^ (^^pl + 12qoVov-\po\^ , (32) 

the width Lo was defined above, and w > is the propa- 
gation velocity. For s = -1-1/ — 1, the first expression rep- 
resents a propagating localized compression/rarefaction, 
while the second denotes a (larger, see comment below; 
cf. Fig. ^ rarefaction/compression, respectively. Notice 
that, for qo ^ a = 0, the first expression recovers the 
KdV result obtained previously (since v jC then recovers 
the KdV soliton width Wi.m), while the second results in 
a divergent (physically unacceptable) solution "^Tl. Fol- 
lowing Wadati, we may re-arrange H3U|I and (|31|l as 



The Extended Korteweg 
Equation 



de Vries (EKdV) 



Repeating the procedure which led to Eq. (|25|) . in 
the previous section, yet now keeping the fourth order 
derivative coefficient q 7^ in Eq. p3|l . one obtains the 
EKdV Equation 



(29) 



where all coefficients are given in ((2611 except a = 
qo/(2wo); recall that a, b are positive by definition. We 
shall see below that po,qo > for Debye interactions (yet 



2V66 
'Sei — X 



9 . 

9^Uan- 



W^^^ tanh 



provided that d ^ qo =/= 0. Here 



W, 



VT - Co 



1/2 



(33) 



(34) 



Furthermore, f; > is the propagation velocity in the 
frame {C,t} and j = 1 (2) recovers ("'2^'') above, so 
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that ei (€2) is equal to +1 (—1), representing rarefactive 
(compressive) solutions, for s = +1 - e.g. the Debye case 
- and vice versa for s = — 1 (which recovers Wadati's 
notation). The pulse width now depends on both Lq 

(defined as previously) and ^'2"'^ The pulse value for 
s — —1 ^3 satisfies: 



= -( + 6av + a) < w^^^ < < w!^^ 



< (Va^ + 6av — a) = w+ , 

(for s = +1, one should permute the superscripts 1 and 
2); since > \w+\, one expects, for s = —1, a small 
rarefactive and a large compressive pulse; the opposite 
holds for s = +1, e.g. in a Debye crystal case: see Fig. 

m 

Inverting to the lattice displacement coordinate u ~ 
J wd(^, expressed in the original coordinates {x,t}, we 
obtain 



{x,t) 



tan 




W^" tanh 



X — Vt — Xq 



where 



{3) 



f \Jv\ + ISgpZiQ {v - vo)-e-i \po\ 
\ \/Po+ 12goWo (w - vq) + ej \po\ 



1/2 



(35) 



(36) 



and j — 1,2. As expected, for any given s {— ±1), the 
two different kink/antikink solutions obtained for differ- 
ent j {— 1 or 2) are not symmetric; cf. Figure |H1 No- 
tice that the maximum value now also depends on W2 
(j = l,2). 

In conclusion, the Extended KdV equation provides a 
more complete description of the nonlinear dynamics of 
the lattice, compared to the KdV equation. In partic- 
ular, the EKdV compressive (rarefactive) pulse soliton 
obtained for s = -1-1 (s = — 1), i.e. po > [po < 0) is 
slightly smaller than its KdV counterpart (see Fig. IHJ, 
but the EKdV also predicts the possibility for a rarefac- 
tive (compressive) soliton, in either case, to form and 
propagate in the same lattice. In the particular case 
of Debye crystals, the net new result to be retained is 
the prediction of the existence of a rarefactive new ex- 
citation, in addition to the rarefactive one, observed in 
experiments. Nevertheless, theoretical studies on molec- 
ular chains seem to suggest that the additional shock-like 
localized mode predicted by the EKdV equation will not 
be as stable as its (KdV-related) counterpart. This pre- 
diction should, therefore, be confirmed numerically (and 
experimentally) before being taken for granted. 



B. The Modified Korteweg - de Vries (MKdV) 
Equation 

Note, for the sake of rigor, that upon setting pq ^ a = 
in Eq.(|25I) above, one obtains a modified KdV (MKdV) 



equation (with only a cubic nonlinearity term). The 
MKdV equation shares all the qualitative properties of 
the KDV Eq. and is, in fact, related to it via a Miura 
transformation (55|. It has two (both negative and pos- 
itive, for each value of s) pulse soliton solutions which 
follow immediately from the preceding solutions H3()|l and 
H31|l of the EKdV equation, upon setting po — 0. The 
remarkable additional aspect of the MKdV equation is 
that it also bears slowly oscillating solutions, named 
breathers, obtained just as rigorously via the inverse scat- 
tering method [2^ l60j| . These solutions (whose wave- 
length is comparable to their localized width, hence the 
'breathing' impression and the name) share the remark- 
able properties of solitons; in particular, they are seen 
to survive collisions between themselves and with pulse 
solitons l22j| . Their analytic form can be readily found 
in Refs. [23| (see §4.1 therein) and ^60] and will not 
be reproduced here, since their condition of existence, 
namely po <C qo (cf. Ref. |22]) is rather not satisfied in 
the case of Debye-inter acting dust grains. Note however 
that breather-like excitations may exist in a DP crystal, 
as one may see via a different (perturbative) analysis of 
the nonlinear modulation of the amplitude of longitudi- 
nal lattice waves. This is considered in separate work 

M. 



VI. THE BOUSSINESQ (BQ) AND 
GENERALIZED BOUSSINESQ (GBQ) 
EQUATIONS 

Remember that the KdV-type equations in the preced- 
ing Section were obtained from the equation of motion 
(|13|l in an approximative manner, assuming near-sonic 
propagation and neglecting high-order time derivatives. 
Those results are therefore expected to hold for veloc- 
ity values only slightly above vq. We shall now see how 
these assumptions can be relaxed by directly relying on 
the initial (nonlinear) equation . 

Let us consider Eq. (|13() again (for = 0). Upon 
setting Po = — 2p, qo = 3q, vfr^ = h > 0, and integrating 
once, with respect to x, one exactly obtains, for w — u^, 
the GENERALIZED BOUSSINESQ (GBq) Equation 



VqWxx ^hWxxxx + P[w^)xx + qi'W^)x 



(37) 



which, neglecting the cubic nonlinearity coefficient q [viz. 
go = in ifT^ ]. reduces to the well-known Boussinesq 
(Bq ) equation, widely studied e.g. in soHd chains; see e.g. 
|2^ I23I Is^. It possesses well-known localized solutions, 
whose derivation is straightforward and need not be re- 
produced here. The exact expressions obtained from l|?7|) 
for the relative displacement w{x, t) and the longitudinal 
displacement u{x,t) are exhaustively presented and dis- 
cussed in Refs. 22, 1231 . The analytic kink/antikink- type 
localized solutions for u{x, t) read 



U3,{x,t) — ^2sgn{h) — I x 

.10. 



6h 



1/2 



10 



tan 



W3 tanh 



vt — Xo 



(38) 



Here the soliton velocity is v, while the soliton width 
depends on both W3 and L3, which are 



[pl + 6qoiv^~vi)f^^T\pa\ 
[pl + 6qo{v^-v'o)Y/^±\po\ 

1/2 



1/2 



(39) 



Recall that, for Debye interactions, h,q > and 
p — —po < (see above), prescribing the 'supersonic' 
(v > vq) propagation of the solutions; the same was true 
of the KdV solitons obtained above. Notice, however, 
that the expressions obtained here for the longitudinal 
displacement u represent both rarefactive and compres- 
sive lattice excitations even for po > (see Table I in , 
for p = —po < 0); remember that this feature was absent 
in the KdV equation (|25|l . where po > i.e. s = +1 al- 
ways led to a compressive solution in H28|l. In fact, this is 
also true of the Bq equation, which is obtained for go = 0, 
i.e. neglecting the last term in the continuum equation 
of motion (|13|l . The exact solution then reads 



UBq{x,t) = -sgnih) sgn{po) 



tanh 



X ~ vt ~ Xq 



(40) 



which, for positive h and poj prescribes only compressive 
supersonic kinks, pretty much like the solution ui derived 
from the KdV theory above (cf. Table I in |23|, for ha > 
0, p < and q = 0). 

Closing this section, one may wish to compare the 
GBq and Bq solutions (|38|l and H40|l , to the homologous 
EKdV- and KdV-related solutions 1(23) and |(2H1) , respec- 
tively, obtained previously: one may readily check that 
the former ones tend to the latter two as v tends to vq 
[to see this, one may set — — {v + 'Vq){v'^ ~ "^o) ~ 
^v^iv^ — uo); recall that h = v\r^. Nevertheless, this 
velocity range restriction is relaxed in the Boussinesq- 
related description. 



VII. EXCITATIONS IN REAL DP CRYSTALS 

It is now quite tempting to observe and compare the 
predictions furnished by the above nonlinear models in 
a DP crystal in terms of excitation features, e.g. dimen- 
sions and form. For instance, one may substitute the 
expressions for the model's physical parameters (i.e. cjqj 
'Vo^ fi, Po and go) derived in illl CI into the definitions in 
the latter three sections, in order to derive a final form 
for localized excitations in a real DP crystal, in terms of 
the propagation velocity w, the lattice parameter k and, 
generally, the sign s (= -1-1 for Debye interactions). The 
interest in this procedure is evident, since one may seek 



feedback (e.g. parameter values, excitation behaviour) 
from experiments and then investigate the validity of the 
above models by adjusting them to real DP crystal val- 
ues. 

The final expressions for Uj{x,t) [j = 1,2,3, cf. H28|l . 
and ||2SJl, respectively] are somewhat lengthy and 
need not be reported here (since they are straightfor- 
ward to derive). We may nevertheless summarize some 
interesting numerical results. 

The soliton width Li, as defined in (|28|l . now becomes 
Li = Vq/poUi^„i- we see that the product of the displace- 
ment u kink's width and maximum value remains con- 
stant (regulated by the cubic interaction potential non- 
linearity), viz. ui,m Li = 1/po, unlike the KdV pulse 
soliton for the relative displacement w — which is 
characterized by i«i,m^i = est.. Both the kink maxi- 
mum value Mi.m and width Li depend on the velocity v; 
as a matter of fact, faster kink excitations will be taller 
and narrower - see Fig. since now 



ro 



Pq 



hi 



1 



v/6(Af-l) 



(41) 



Recall that the Mach number M — v/vo is always larger 
than unity. Furthermore, the magnitude of the excitation 
seems to decrease with k; see Fig|^: nevertheless, very 
high values (near u/ro = 1) observed for low k and high 
V are rather not to be trusted, since they contradict the 
continuum approximation tt <C tq. 

Finally, one may compare the solutions obtained from 
the above theories, for a typical value of k, say 1.25, 
according with real experimental values. The KdV-, 
EKdV- and Boussinesq-related kink excitations, i.e. mi, 
U2 and U3, are depicted in Figs. ^1-^1 for three differ- 
ent values of M, 1.1, 1.25 and 2. We see that the EKdV 
and Bq models allow for both compressive and rarefac- 
tive structures, while the KdV description predicts a lo- 
calized compression, which is quite sensitive to velocity 
changes. As expected (cf. the discussion above), both 
the EKdV- and Bq-related compressive excitations are 
similar in magnitude to the KdV-related anti-kink for 
near-sonic velocity (i.e. near M « 1). Nevertheless, we 
see that the KdV-related antikink becomes taller and 
narrower as velocity increases, and substantially differ- 
entiates itself from its EKdV- and Bq-analogues. One 
may wonder whether or not the KdV picture (more fa- 
miliar since widely studied) is adequate for the modeling 
of a real DP crystal, and also whether the rarefactive ex- 
citations predicted by other theories can indeed be sus- 
tained in the crystal. These questions may be answered 
by appropriate experiments and, possibly, also be inves- 
tigated by numerical simulations. From a purely theo- 
retical point of view, the Boussinesq-based description 
appears to be more rigorous (recall that the KdV was 
derived in some approximation) and valid in a more ex- 
tended region than both the KdV and Extended-KdV 
theories. 



11 



VIII. DISCRETENESS EFFECTS 

The above analytical solutions have been derived in 
the continuum limit, i.e. for L » tq, where L is the 
typical spatial dimension (width) of the solitary excita- 
tion. One may therefore define the discreteness parame- 
ter g = r^l L, and require a posteriori that g From 
the expressions derived for the Bq equation above, one 
easily sees that g ~ (u^ — Uq)^/^/wi, so this requirement 
is indeed fulfilled for propagation velocities u « . How- 
ever, for higher values of v, the (narrower) soliton will 
be subject to a variety of effects e.g. shape distortion, 
wave radiation etc., due to the intrinsic lattice discrete- 
ness. These effects have been investigated in solid state 
physics and may be considered with respect to 

DP crystals at a later stage. Let us briefly point out that 
narrow kink-shaped lattice excitations have been numer- 
ically shown to propagate with no considerable loss of 
energy, in a quite general monoatomic lattice model [62| . 

Also worth mentioning is the work of Rosenau [g^ who 
derived an improved version of the Boussinesq equation 
(the I-Bq Eq.) in a quasi-continuum limit. The I-Bq 
equation, which bears the general structure of l|37|) upon 
replacing h u^xxx therein by h Uxxtt (yet with different co- 
efficient definitions), is not integrable and bears solitary 
wave solutions which do not collide elastically; neverthe- 
less, it was numerically shown to be more stable than the 
Bq equation, and was argued to model discrete lattice dy- 
namics more efficiently, upon comparison of theoretical 
predictions to exact numerical results . Further exam- 
ination of such effects may be carried out in dust-lattices, 
once our feedback from experiments has sufficiently de- 
termined the relevance of the issue in real DP crystals, 
i.e. typical excitation width, dynamics etc. 

It should be underlined that the possibility for the exis- 
tence of breather solitons, anticipated above, establishes 
a link between complex plasma 'solid state' modeling and 
the framework of discreteness-related localized excita- 
tions (discrete breathers [g^, intrinsic localized modes 
|67|). which have recently received increasing interest 
among researchers in the nonlinear dynamics community. 
These localized modes, which are due to coupling anhar- 
monicity and are stabilized by lattice discreteness, have 
been shown to exist in frequency regions forbidden to or- 
dinary lattice waves and account for energy localization 
in highly discrete real crystals, where continuum theories 
fail. The relevance of this framework to dust crystals ap- 
pears to be an interesting open area for investigation. 



IX. ENVELOPE EXCITATIONS AND SHOCKS - 
OPEN ISSUES 

As a final interesting issue involved in the nonlinear 
dynamics of longitudinal lattice oscillations, let us men- 
tion the nonlinear modulation of the amplitude of dust- 
lattice waves, a well-known mechanism related to har- 
monic generation and, possibly, the modulational insta- 



bility of waves propagating in lattices, eventually lead- 
ing to modulated wave packet energy localization via 
the formation of envelope solitons [23 ■ This framework, 
which was recently also investigated with respect to low- 
frequency (dust-acoustic, dust-ion acoustic) electrostatic 
waves in dusty plasmas [s^, has been partly analyzed, 
on the basis of the Melands0[l5| model in Ref. Is^]. The 
authors relied on a truncated Boussinesq equation, in 
the form of H37() for q — 0, and succeeded in predict- 
ing the occurrence of modulational instability in LDL 
waves in DP crystals and the formation of envelope struc- 
tures. Nevertheless, the nonlinearity coefficient q omitted 
therein seems to compete with p in H37() (notice the dif- 
ferent signs) and is rather expected to affect significantly 
the wave's stability profile. It should be stressed that 
these localized envelope excitations result from a phys- 
ical mechanism which is intrinsically different from the 
one related to the small-amplitude excitations described 
in this paper; see the discussion in Ref. 68]. An ex- 
tended study of this modulation nonlinear mechanism is 
in row and will be reported elsewhere, for clarity and 
conciseness. 

As a final comment, we may speculate on the role of 
damping, herewith ignored, on the dynamics of dust- 
lattice waves. It is known that weak damping may bal- 
ance nonlinearity, leading to the formation of shock wave 
fronts, as predicted in Refs. .42J and confirmed by 
numerical simulations Furthermore, it was recently 
shown that the same mechanism may result in the for- 
mation of large-amplitude wide-shaped solitary waves, 
which may later break into a (gradually damped) train 
of solitons or a wavepacket depending on physical pa- 
rameters We see that friction, yet weak, may play 
a predominant role in the life and death of localized ex- 
citations; this effect definitely deserves paying close at- 
tention with respect to waves propagating in dust crys- 
tals. Again, one would expect phenomenological theories 
followed by appropriately designed experiments to elu- 
cidate the friction mechanisms inherent in longitudinal 
dust-lattice wave propagation, in view of a more com- 
plete description than the one provided by the conserva- 
tive model adopted here. 



X. CONCLUSIONS 

This work was devoted to an investigation of the rel- 
evance of existing model nonlinear theories to the dy- 
namics of longitudinal oscillations in anharmonic chains, 
with emphasis on dust-lattice excitations in (strongly- 
coupled) complex plasma crystals. Taking into account 
an arbitrary interaction potential and long-range inter- 
actions, we have rigorously shown that both compres- 
sive and rarefactive kink-shaped (shock-like) excitations 
may form and propagate in the lattice, depending basi- 
cally on the mechanism of interaction between grains lo- 
cated at each site. These excitations are effectively mod- 
eled by either KdV- or Boussinesq-type equations, whose 



12 



analytic form was presented and whose qualitative and 
quantitative differences were discussed. In any case, the 
theory predicts coherent wave propagation above the lat- 
tice's 'sound' speed, in agreement with previous theoreti- 
cal works and experimental observations (in both atomic 
and dust- lattices). It may be appropriate to mention 
that subsonic soliton propagation in monoatomic chains 
was also numerically considered and shown to be feasible 
in the past jo^. Let us point out that the model used 
here to pass from a discrete description to the continuum 
(long-wavelength) limit is quite generic, so possible mod- 
ification via refined nonlinear equations may readily be 
obtained from it, for future consideration. 

Furthermore, we have discussed the possibility of the 
formation of breather modes and envelope excitations, as 
a consequence of modulated wave packet instability, an- 
ticipating their link to discrete nonlinear theories of lo- 
calized modes, left for future consideration; despite their 
analytical complexity, these models may, in principle, be 
of relevance in dust crystals due to the finite dimensions 
of the chain and its intrinsic spatial discreteness. Fi- 
nally, the possible role played by dissipation mechanisms 



has been briefly discussed. 

The present study relies on, and aims at extending, 
previous theories on both anharmonic atomic chains and 
dusty plasma crystals. We hope to have succeeded in re- 
viewing the former (extending them to the case of long- 
range electrostatic interactions) and generalizing the lat- 
ter (which are still in an early stage) . Hopefully, our pre- 
dictions may be confirmed by appropriately set-up exper- 
iments, with the ambition of throwing some light in the 
relatively new and challenging field of strongly- coupled 
complex plasmas and dust-lattice dynamics. 
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APPENDIX A: COMPUTATION OF THE 
COEFFICIENTS FOR DEBYE (YUKAWA) 
INTERACTIONS 

Consider the Debye potential (pDif) = Qe^^'l^^jr. 
Defining the (positive real) lattice parameter k = r^/X^, 
it is straightforward to evaluate the quantities 



I + Ik, 



2g + + ^ 



{In) 



24Q i + + 



24 



(In 



A5 



where the prime denotes differentiation and I = 1, 2, 3, ... 
is a positive integer. Now, we shall combine these expres- 
sions with Eqs. iQ, ||SJ), and defining Vg, vf, 
Po and qo, respectively. 

Let us define the general (families of) sum(s) 

oo N 

5„ (a) = ^ a' r Si^^ (a) = ^ a' ^" 

1=1 1=1 

(0 < a < 1) , (Al) 

(thinking of a = e"", in particular); note that Sn^\a) — > 

Sn{a) for ^ oo; also, Sn\a) = a. Making use of the 
well-known geometrical series properties: 



N 



1=1 



1=1 



then 



^3 (a) 



^3 g2 Q 



a(a2 + 4a + 1) 
^ - (l-a)4 

and so forth. Also note the identity 

J 



1=1 



S-i{a) = E] — = —ln{\ — a) 



(0 < a < 1) . 
(A3) 

The corresponding set of formulae may be obtained for 
Sri^' in a similar manner. 

Now, substituting a = e"" and using the derivatives of 
above, one may immediately evaluate the expressions 
O, (IHl, lUni) and (inji. setting ro = nXo everywhere, it 
is straightforward to show that 

— 2_2 2 Q 2 \2 i2 illn \ 

C2 = Vq = Wq.l^o = ]^ ^dZ^I <f> [Iro) = ... 



1=1 



9(1'^ r 1 



(A4) 



In the same manner 



C4 _ 2 



Q 

12M 



K'\lY.l'c^"{lr,) 



1=1 



GMXr 



I _ a(l-a^) Also 



(0 < a < 1) , (A2) 

it is straightforward to derive Sn, Si^'^ for /> 1, by 
differentiating with respect to a. One obtains 

oo oo oo /X 



Po 



K-^Si{e~'') + K°52(e-") 

+ ^^^'^3(6"=) 



^ -cii =-%^'XlY.l'^"'(l'o) 



(A5) 



M 



1=1 



6Q2 



1=1 



d 



1=1 



1=1 



do 



, dSo a 
a — a- 



(A6) 



" 5a ^ da (1 - ay 

In a similar manner, iterating from 

d^{a^)/da'^ = a'-2 = (;2^/ _ 

one finds 



Finally, from p2|) we have 

Q 



2M 



n''X%J2l'r'(.lro) = 



1=1 



a(l -I- a) 
(T^ 



go = Cm 

^ 12g^ 
~ MXd 



K-^S-iie-'') + K°^o(e'") 



(A7) 
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The corresponding expressions for a value of N are given 
by substituting £'„(•) with S^\-) everywhere. One im- 
mediately sees that po/vq > 2, qo/vq > 6; also, vf/vQ = 
12 for iV = 1 (only), i.e. for the NNI case. 

Finally, combining the above exact expressions for 
S-i{a), 5*3(0), we obtain exactly expressions (fl^ to 
in the text. 



Figure Captions 

Figure 1. 

(a) The linear oscillation frequency squared cj^ (nor- 
malized over Q^/(Af A|,)) is depicted against the lattice 
constant k, for = 1 (first-neighbor interactions: — 

), N = 2 (second- neighbor interactions: ), N = 00 

(infinite- neighbors: ), from bottom to top. (b) De- 
tail near k « 1. 

Figure 2. 

(a) The characteristic 2nd order dispersion velocity 
squared Vq (normahzed over Q^/(MA_d)) is depicted 
against the lattice constant k, for = 1 (first-neighbor 
interactions: — ), N = 2 (second-neighbor interactions: - 

- -), iV = cxD (infinite-neighbors: ), from bottom to 

top. (b) Detail near k « 1. 

Figure 3. 

(a) The characteristic 4th order dispersion velocity 
squared vf (normahzed over Q'^ /{M\d)) is depicted 
against the lattice constant k, for = 1 (first-neighbor 
interactions: — ), N = 2 (second-neighbor interactions: - 

- -), N ~ 00 (infinite-neighbors: ), from bottom to 

top. (b) Detail near k w 1. 

Figure 4. 

(a) The nonlinearity coefficient po (normalized over 
(5^/(MAd)) is depicted against the lattice constant n for 

= 1 (first- neighbor interactions: — ), N = 2 (second- 
neighbor interactions: ), N = 00 (infinite-neighbors: 

), from bottom to top. (b) Detail near k sa 1. 

Figure 5. 

(a) The nonlinearity coefficient qq (normalized over 
Q"^ /{M\d)) is depicted against the lattice constant k for 
A^ = 1 (first-neighbor interactions: — ), N = 2 (second- 
neighbor interactions: ), N = 00 (infinite-neighbors: 

), from bottom to top. (b) Detail near k « 1. 

Figure 6. 

Dispersion relation for the Debye interactions, neglect- 
ing damping; cf. for v — 0: the square frequency 
Lu'^, normalized over Q'^/{MXy), is depicted versus the 
normalized wavenumber kro/ir for A'^ = 1 (first-neighbor 
interactions: — ), N — 2 (second-neighbor interactions: 

), N — 7 (up to 7th nearest-neighbors: ), i.e. 

from bottom to top. 

Figure 7. 

Localized antikink/kink (negative/positive pulse) 
functions, related to the KdV Eq. H25I) . for the 
displacement u{x,t) (relative displacement 'w{x,t) ~ 
du{x,t)/dx), for positive/negative po coefficient i.e. s = 
-1-1/ — 1, are depicted in figures (a) /(b); recall that (a) 
holds for Debye interactions; arbitrary parameter values: 
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V = 1 (solid curve), 
(short dashed curve). 



2 (long dashed curve), v = 3 



Figure 8. 

(a) The two locaUzed pulse solutions of the EKdV Eq. 
(|29|l for the relative displacement 'w{x,t) ~ du{x,t) /dx 
are depicted for some set of (positive) values of the 
po and qo coefficients (i.e. s = +1): the first 
(dashed curve)/second (short-dashed) solution, as given 
by (|30|l / H31() . represents the smaller negative/larger pos- 
itive pulses. The larger negative pulse (solid curve) de- 
notes the solution of the KdV Eq. H25I) for the same 
parameter set. (b) The corresponding solutions for the 
particle displacement u{x,t). 
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Figure 9. 

(a) The (normalized) maximum value of the kink- 
shaped localized displacement ui{x,t)/ro, as obtained 
from the KdV Equation, is depicted versus the lattice 
parameter k and the normalized velocity (Mach number) 
M = v/vq. (b) The (normalized) width Li/tq of ui{x,t) 
is depicted against M = v/vq. 

Figure 10. 

The antikink excitation predicted by the KdV theory 
(solid curve) is compared to the (two) solutions obtained 
from (a) the EkDV Equation; (b) the Bq Equation- 
related model (dashed curves). Values: lattice param- 
eter K = 1.1, normalized velocity (Mach number) M = 
v/vq = 1.25. 
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Figure 11. 

Similar to Fig. [TOI for M = v/vq = 1.25. 



Figure 12. 

Similar to Figs. [n]and[T21 for M = v/vq = 2. 
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sound velocity vO''2 vs k 
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nonlinearity coeff. pO vs kappa nonlinearity coeff. qO vs kappa 
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dispersion relation omega'^2 vs k 
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KdV-related solutions , s = +1 
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KdV vs . Boussinesq, M = 2 
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